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ABSTRACT 

We examine the amplitude of the rotation velocity that can be attributed to 
the dark matter halos of disk galaxies, focusing on well measured intermediate 
radii. The data for 60 galaxies spanning a large range of mass and Hubble types, 
taken together, are consistent with a dark halo velocity log Vh = £ + 23 log r 
with £ = 1.47^o!i9 an d 23 ~ \ over the range 1 < r < 74 kpc. The range in 
£ stems from different choices of the stellar mass estimator, from minimum to 
maximum disk. For all plausible choices of stellar mass, the implied densities of 
the dark halos are lower than expected from structure formation simulations in 
ACDM, which anticipate € > 1.6. This problem is not specific to a particular 
type of galaxy or to the innermost region of the halo (cusp or core); the velocity 
attributable to dark matter is too low at all radii. 

Subject headings: dark matter — galaxies: kinematics and dynamics — galaxies: 
spiral 



1. Introduction 

The rotation curves of spiral galaxies provide one of the strongest lines of evidence 
establishing the need for dark matter in the universe (e.g., Rubin, Thonnard & Ford 1980; 
Bosnia 1981). An elaborate paradigm for dark matter has subsequently developed: ACDM. 
More recent work emphasizes larger scale constraints (e.g., Tegmark et al. 2004; Spergel et 
al. 2006). An important question is whether this cosmic paradigm provides a satisfactory 
description of spiral galaxies. 
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A complete explanation of spiral galaxies in ACDM requires a comprehensive theory of 
galaxy formation. This remains lacking. Indeed, there are a number of lingering problems 
on small (galaxy) scales. These include the cusp-core problem (e.g., de Blok et al. 2001a), 
the missing satellite problem (e.g., Moore et al. 1999a), the dynamical friction time scale 
problem (Goerdt et al. 2006; Sanchez-Salcedo, Reyes-Iturbide, & Hernandez 2006), and a 
whole suite of other incongruities and apparent contradictions that arise on the scale of 
individual galaxies (McGaugh & de Blok 1998; Sellwood & Kosowsky 2001). In this paper, 
we will concentrate on the amplitude of the rotation velocity that is observed for spiral 
galaxy dark matter halos and that predicted by ACDM. 

Considerable attention has focused on the density profile inferred for the dark matter 
halos of dwarf and low surface brightness (LSB) galaxies at small radii (Moore 1994; Flores & 
Primack 1994; Cote et al. 2000; Swaters et al. 2000; de Blok et al. 2001a,b; Borriello & Salucci 
2001; Salucci 2001; Bolatto et al. 2002; Marchesini et al. 2002; Salucci et al. 2003; Swaters 
et al. 2003a,b; de Blok et al. 2003; Donato et al. 2004; Simon et al. 2003, 2005; Gentile et al. 
2004, 2005; Spekkens, Giovanelli, & Haynes 2005; Kuzio de Naray et al. 2006; Zackrisson et 
al. 2006). There are good reasons for this: LSB galaxies are dark matter dominated down to 
small radii, making the inference of their dark matter properties much less sensitive to the 
stellar mass than in higher surface brightness spirals. In the most straightforward analyses, 
the data appear to contradict a basic prediction of the ACDM structure formation scenario, 
that dark matter halos should have cusps at small radii (a density distribution p ~ r~ 7 with 
7 w 1: Dubinski 1994; Cole & Lacey 1996; Navarro, Frenk, & White 1997 (hereafter NFW); 
Tormen, Bouchet, & White 1997; Moore et al. 1998, 1999b; Tissera & Dominguez-Tenreiro 
1998; Nusser & Sheth 1998; Syer & White 1998; Avila-Reese, Firmiani, & Hernandez 1998; 
Salvador-Sole, Solanes, & Manrique 1998; Jing 2000; Jing & Suto 2000; Kull 1999; Klypin 
et al. 2001; Ricotti 2003; Power et al. 2003; Diemand et al. 2005). The data are apparently 
more consistent with a nearly constant density core (7 w 0). Many different solutions to this 
problem have been proposed (e.g., El-Zant, Shlosman, & Hoffman 2001; Weinberg & Katz 
2002; Merritt et al. 2004; Ma & Boylan-Kolchin 2004; Mashchenko, Couchman, & Wadsley 
2006; Tonini, Lapi, & Salucci 2006). 

Interesting constraints also stem from the rotation curve data for high surface brightness 
(HSB) galaxies. It is frequently possible to obtain quite high quality data for these bright 
objects, but the interpretation depends critically on the mass ascribed to baryons. The 
stellar disks in these objects can be massive (e.g., Debatista & Sellwood 2000; Salucci & 
Burkert 2000; Weiner, Sellwood, & Williams 2001; Zavala et al. 2003; Kassin, de Jong, & 
Weiner 2006) or even maximal (e.g., van Albada & Sancisi 1986; Kent 1987; Palunas & 
Williams 2000; Battaglia et al. 2006). Such dynamically dominant disks leave little room 
for halos with steep inner density profiles. For example, Binney & Evans (2001) find that 
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cuspy halos are inconsistent with a number of constraints for the Milky Way (cf. Klypin, 
Zhao, & Somerville 2002). On the other hand, it can be argued that the lack of scale length 
residuals in the Tully-Fisher relation implies that typically HSB disks are sub-maximal (e.g., 
Courteau & Rix 1999; cf. Sellwood 1999). If dark matter is dynamically significant in the 
inner regions of HSB galaxies, then it may be quite possible for the halos of these galaxies 
to have a steep cusp. The critical issue is the stellar mass. 

Constraining the precise inner slope of the dark matter density profile of an individual 
galaxy is observationally challenging. Even if the stellar mass can be properly treated, a 
variety of systematic effects might affect the result (Swaters et al. 2003a; de Blok et al. 
2003; Simon et al. 2003, 2005; Spekkens et al. 2005). There are several distinct problems 
here: observational uncertainties, non-circular motions that might cause misestimation of 
the potential, and degeneracy between fitting parameters. 

Observational uncertainties of concern here are those that might lead us to underesti- 
mate the steepness of the rise of the rotation curve. This might result, for example, from 
beam-smearing of 21 cm observations (e.g., van den Bosch et al. 2000) or mis-positioning of 
the slit during optical observations (e.g., Swaters et al. 2003a). Such errors are not a great 
concern when rotation curves are observed to be steeply rising, as is typically the case for 
HSB galaxies in this sample. When the observed rotation curve rises more gradually, there is 
more potential for concern as to whether this behavior is truly intrinsic (de Blok, McGaugh, 
& van der Hulst 1997; de Blok & McGaugh 1997). This no longer appears to be a significant 
worry for LSB galaxies, as high resolution two dimensional velocity fields (Kuzio de Naray 
et al. 2006) largely confirm pervious results (e.g., McGaugh, Rubin, & de Blok 2001; de Blok 
& Bosma 2002; see also Zackrisson et al. 2006). 

Non-circular motions can be significant at small radii, in the sense that they can repre- 
sent a substantial fraction of the total gravitational potential. The contributions of circular 
and noncircular motions sum in quadrature. From this, it immediately follows that non- 
circular motions (V r ) are only important when comparable to the rotational motion (V^). 
Simon et al. (2005) show some cases for which this is the case (V r « V^), albeit at quite 
small radii. For the sample considered here, which all have two dimensional velocity fields, 
V r is typically small (a few to tens of kms" 1 ) compared to at the same radius (tens to 
hundreds of kms -1 ). We consider only intermediate and large radii where S> and the 
contribution of non-circular motions to the potential is not significant. Galaxies that do not 
satisfy this criterion are not in the sample. 

A more fundamental difficulty is the degeneracy between halo fitting parameters, even 
when the potential is perfectly well known. The degeneracy between stellar mass and the core 
radius of a pseudo-isothermal halo has long been known (e.g., Kent 1987). For NFW halos, 
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there is substantial covariance between c and V200 even when the stellar mass is fixed (de 
Blok et al. 2001b). This degeneracy stems from the similarity of the inner shapes of NFW 
halo velocity profiles over a range of (c,V 2 oo)- More generally, degeneracy is unavoidable 
because of a surplus of parameters. Mass models require the specification of a minimum of 
three parameters per galaxy (one for the disk and two for the halo) while empirically a single 
parameter suffices to describe the data (e.g., Brada & Milgrom 1999; McGaugh 2004). 

In this paper, we adopt a model independent approach. We do not fit specific mass 
models to individual galaxies. The halo parameters derived from such models inevitably 
suffer from parameter degeneracies. Rather, we simply ask what range of halo parameters 
(c, V200) are consistent with the halo velocities implied by observation. We determine the 
velocity attributable to the dark matter for a broad range of stellar mass estimators. 

We find that the data are consistent with a nearly universal dark matter halo profile. 
There is some scatter about this relation, and we would not claim that there is in fact 
a single halo profile appropriate to all galaxies, but the data are not wildly at odds with 
such a situation. The slope we find is consistent with previous determinations (e.g., Barnes, 
Sellwood, & Kosowsky 2004). If halos have the NFW form, their concentrations and/or 
masses (V 2 oo) are generally smaller than might nominally be expected in ACDM. This is 
largely consistent with various previous results (e.g., Debatista & Sellwood 2000; Weiner et 
al. 2001; Borriello & Salucci 2001; Zentner & Bullock 2002; Alam, Bullock, & Weinberg 2002; 
McGaugh et al. 2003; Bell et al. 2003a; Gentile et al. 2004; Zackrison et al. 2006; Portinari & 
Sommer-Larsen 2006) that suggest that real dark matter halos are less dense than predicted 
by ACDM. 

The data and stellar mass estimators are described in §2. In §3 we determine the halo 
velocity V/, implied by the data for a range of stellar mass estimates. The predictions of 
ACDM models and their relation to the data are described in §4. The implications of the 
result and possible ways to reconcile the apparent discord between theory and observation 
are discussed in §5. A summary of our conclusions is given in §6. 

2. Data 
2.1. Sample 

We employ the rotation curve data and mass models compiled by Sanders & McGaugh 
(2002), as trimmed for quality by McGaugh (2005). These data originate from many in- 
dependent observers, and represent the accumulated effort of a substantial community of 
astronomers. Sample galaxies span a wide range of morphological types, from Sab through 
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all spiral types to Irr/Im. Recent work by Noordermeer et al. (2005) focussing on the earli- 
est (SO/Sa) type disk galaxies appears contiguous with the data used here. More important 
than morphological type is the broad range of physical parameter space that our sample 
spans. It includes galaxies with circular velocities 50 < Vf < 300 kms -1 , baryonic masses 
3 x 10 8 < Aib < 4 x 10 11 M.Q., scale sizes 0.5 < Rd < 13 kpc, central surface brightnesses 
19.6 < < 24.2 mag. arcsec -2 , and gas fractions 0.07 < f g < 0.95. No restriction is made 
on surface brightness or any other physical parameter. The sample represents the dynamics 
of disks in the broadest sense possible with available data. 

All sample galaxies have high quality, extended rotation curves derived from 21 cm 
data cubes. The high precision sample of McGaugh (2005) consists of 60 galaxies with 
~ 700 independent, resolved rotation velocity measurements, each with formal accuracy of 
5% or better. We choose to exclude data with r < 1 kpc where the controversy over cusp or 
core occurs (de Blok 2004). The point here is to compare reliable data with the expectations 
of ACDM models, free from concern over systematic uncertainties, either observational (such 
as beam smearing) or physical (e.g., non-circular motions). These may be serious concerns 
for constraining the innermost slope of the density profile, but do not significantly affect the 
data at intermediate radii. Here the data are well resolved, and the observed circular velocity 
greatly exceeds the velocity dispersion. Indeed, strong constraints exist on the ellipticities of 
the potentials for some sample galaxies (Schoenmakers, Franx, & de Zeeuw 1997). Triaxiality 
of the halo (e.g., Hyashi & Navarro 2006) is not relevant to the bulk of the data here. 

We exclude the possibility of non-circular motions biasing the results at small radii by 
truncating the inner data. Our result is not sensitive to the choice of the inner truncation 
radius. The total radial range sampled is from 1 to 74 kpc, with a median radius of 9 kpc. 
Only a few galaxies have measured rotation extending to very large radii. The bulk of the 
data (68%) fall in the range 4 < r < 22 kpc (Fig. [T]). Data in this range are not only 
numerous, but fall in a range where the systematic concerns discussed above have negligible 
impact. The same result follows if we restrict the data to this very conservative range. 

For many galaxies in the sample, the contribution of the stars to the total potential 
can not be ignored. This dominates all other effects by far. It is necessary to subtract off 
the baryonic contribution to reveal the velocity due to the halo. All galaxies have detailed 
surface photometry from which baryonic mass models are derived. 

For the present work, we utilize -B-band surface photometry, as this is the only band 
that all 60 galaxies have in common. We are working to place the sample on a uniform K s - 
band scale, but at present lack sufficiently deep observations of some of the fainter galaxies. 
Using the i^-band data already in hand does not alter the basic result (see also Verheijen 
2001; Verheijen & Sancisi 2001; Kassin et al. 2006). 
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We make no assumption about the light profile (e.g., exponential disk). The observed 
light distribution is converted to a baryonic rotation curve numerically (via GIPSY: van der 
Hulst et al. 1992). This includes the bulge (when present), stellar disk, and Hi gas. The 
molecular gas distribution is not available in most cases, but this is likely to be a small 
percentage of the baryonic mass in the majority of late type galaxies (Young & Knezek 
1989). Moreover, the molecular gas generally follows the stellar distribution (Heifer et al. 
2003), so is effectively subsumed into the stellar mass-to-light ratio. It is the estimation of 
stellar mass that is most crucial. 



2.2. Stellar Mass 

The observed rotation curves provide a good map of the total gravitational potential 
that is due to both baryonic and dark mass. While the total potential is well constrained, 
dividing it into baryonic and dark mass components is a degenerate process. Additional 
information is needed to break this degeneracy. 

We explore three methods for estimating the stellar mass-to-light ratio Y*: maximum 
disk, T max ; stellar population synthesis models, T pop ; and that from the mass discrepancy- 
acceleration relation, T acc . This latter in effect treats MOND (Milgrom 1983) fits (Sanders 
& McGaugh 2002) as an empirical phenomenology. That this phenomenology is present in 
the data is all that is required, regardless of the underlying physical reason for it (McGaugh 
2004). 

Following McGaugh (2005), we can choose any arbitrary stellar mass estimate by scaling 
each of these prescriptions by a constant factor T, V, or Q: 

y '* rT " max (i^) 
T* = VT pop (lb) 
T* = QT acc . (lc) 

This scheme allows us to explore the full range of possible stellar masses (0 < T < 1) with 
three independent estimators for each galaxy. The population synthesis mass-to-light ratio 
T pop is estimated from the observed B — V color (as tabulated by McGaugh 2005) and the 
models of Bell et al. (2003b): 

logT pop = 1.737(5 -V)- 0.942 (2) 

Very similar results follow from other models (e.g., Portinari et al. 2004). The chief un- 
certainty of this method is the IMF. V — 1 corresponds to the nominal "scaled Salpeter" 
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IMF adopted by Bell et al. (2003b), which for our purposes is nearly indistinguishable from 
the Kroupa (2002) IMF adopted by Portinari et al. (2004; see also Portinari & Sommer- 
Larsen 2006). To consider the effects of a heavier or lighter IMF, we simply scale by V (see 
discussion in Bell & de Jong 2001). 

3. Empirical Halo Estimates 

We estimate the circular velocity provided by dark matter by subtracting the baryonic 
component from the observed rotation curves for each choice of Y*: 

V^r) = V\r)-V h \r). (3) 

Here V^(r) is the rotation velocity attributed to the dark matter, V(r) is the observed 
rotation velocity, and V&(r) is the rotation velocity attributed to baryons. This depends on 
Y* through V b 2 = T+v^ + V g 2 , where v± is the rotation velocity for Y* = 1 A4q/L q . The 
Hi gas component, corrected to account for helium and metals, gives V g . Molecular gas 
generally traces the stellar distribution (Heifer et al. 2003) and is subsumed into T*. 

An example rotation curve is shown in Fig. [2j As we increase (decrease) Y* by varying 
T, V, and Q, the stellar contribution to the total rotation increases (decreases) and the 
corresponding contribution from the dark matter decreases (increases). Though the details 
of each realization depend on the T* estimator and the relative importance of gas (which 
does not vary), it is a tolerable approximation to imagine the dark and baryonic lines in 
Fig. [2] sliding up and down see-saw fashion. 

In order to restrict our attention to the best data, we have excised individual points 
with uncertainty in excess of 5%. In many cases, the formal precision is rather greater than 
this limit. While this leaves plenty of data to analyze, various astrophysical considerations 
make it unrealistic to consider those data with small error bars as accurate estimates of the 
circular velocity of the potential (see discussions in McGaugh et al. 2001; de Blok et al. 
2003). We therefore give uniform weight to the data which pass the 5% criterion. 

A more serious uncertainty than that in the data is that in T*. The estimate of stellar 
mass is the single dominant source of uncertainty in this analysis. The systematic nature of 
this uncertainty makes it impossible to propagate it in a meaningful way. Hence we simply 
examine the results of adopting each Y* estimator, leaving it to the reader's judgement which 
case might seem best. Fortunately, the basic result is fairly insensitive to this choice. 

In Fig. [3] we show Vh (r) for all 60 galaxies plotted together for the case of Q = 1. 
Remarkably, the inferred dark matter halos of all these galaxies fall more or less on top of 
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each other. This is not to say that the dark matter halos of all galaxies are the same. Quite 
the contrary, it is possible to observationally distinguish between the dark matter halos of 
individual galaxies (Sancisi 2004; McGaugh 2004; Barnes et al. 2004; Simon et al. 2005). 
Nonetheless, when plotted together in this fashion, the apparent variation is remarkably 
small. 

Taken together, the data suggest a relation 

logVfc = e+»logr. (4) 

Fitting the data with Q = 1 gives € = 1.48 ± 0.01 and 03 = 0.49 ± 0.01. While we 
would not claim that this is a universal halo that all galaxies share, it does give a tolerable 
approximation of the velocity attributable to dark matter. 

We can extend this approach to other mass-to-light values by varying T, V, and Q. The 
cases r = 1, V — 1, Q — 1 and T = 0.4 (the latter as suggested^ by Bottema 1993 and 
Courteau & Rix 1999) are illustrated in Fig. HI This figure is like Fig. for the various mass 
estimators, but with shading to illustrate the density of data. The darker the shading, the 
more sharply defined the relation between Vh and radius. 

There is some covariance between the slope 03 and intercept <£. Our purpose here is 
to constrain the amplitude of the velocity that can be attributed to dark matter. In order 
not to obscure the basic result with uncertainties introduced by covariance, we fix the slope 
to 03 = 1/2, consistent with our data (see Figs. [5] -[8]) and those of Barnes et al. (2004). 
Moreover, the slope does not vary appreciably when left free; the important information is 
contained in the amplitude. 

Fixing 03 = 1/2 changes the Q = 1 best fit intercept to € = 1.47. This number encodes, 
in the simplest possible way, the amplitude of the dark matter velocity. Fig. [9] show how 
this amplitude varies as we change T*. For V = 1, the amplitude barely changes: £ = 1.49. 
For the full range of possibilities, from minimum to maximum disk, we find €(T = 0) = 1.62 
and £(r = 1) = 1.28. In other words, the amplitude of the velocity provided by the dark 
matter falls in the range £ = 1.47l^g for 03 = 1/2. 

One can also see in Fig. [9] variation in the scatter about the mean. Remarkably, the 
scatter is minimized for Q — 1, the MOND mass-to- light ratios. This is also the choice 
of which minimizes the scatter in the baryonic Tully-Fisher relation (McGaugh 2005) 
and, by construction, that in the mass discrepancy-acceleration relation (McGaugh 2004). 



1 These authors phrase the disk contribution in terms of the velocity at 2.2 scale lengths. Since Ai oc V 2 , 
V2.2 /Vtot = 0.63 corresponds to T = 0.4. 
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Irrespective of the physics that causes it, this implies that MOND is an important empirical 
organizing relation that is relevant even to dark matter. The scatter is larger for other mass 
estimators, but the basic result encapsulated by £ is the same. 



Cosmological simulations lead to mass distributions for dark matter halos that have 
gravitational potentials which can be analytically approximated as 



where c = R.2oo/Rs an d x = r/i? 2 oo (NFW). The radius -R200 encloses a density 200 times 
the critical density of the universe; V200 is the circular velocity of the potential at -R200 

The radius -R200 is, crudely speaking, the virial radius in a critical density universe. For 
current ACDM parameters, the virial radius is more nearly at R\oo (e.g., Eke, Navarro, & 
Steinmetz 2001). More recently, Navarro et al. (2004) have suggested that a density profile 
with a continuously running slope provides a better fit to simulated halos than the original 
NFW formula. From an observational perspective, these forms are virtually indistinguish- 
able. We therefore retain the NFW formulation for its familiarity and its well quantified 
relation to cosmological parameters (NFW; McGaugh, Barker, & de Blok 2003). 

One interesting aspect of the data is that there is no clear tendency for the dark halo 
rotation curves to flatten out. The combination of disk plus halo is often close to flat, and 
as the disk contribution is declining, it is natural to associate the flat rotation velocity with 
the halo. This is, however, an oversimplification. In most cases, the inferred dark matter 
contribution to the velocity is still rising at the last measured point, and it is not clear that 
the observed flat velocity is truly representative of the halo. 

In the context of simulated halos, the lack of a flat halo rotation velocity may help the 
models. NFW halos do not possess a significant range of radii where p cx r~ 2 , as needed to 
produce flat rotation curves. Presumably, the observed range suffices only to cover the rising 
portion of the halo, with a suggestion of a turn-over in only a few sufficiently extended cases. 
It remains an open question why the disk and halo conspire to produce nearly flat rotation 
curves over the range observed. 



4. 



ACDM Halos 




ln(l + cx) — cx/(l + cx) 



(5) 



It is suggestive that the fitted slope OS ~ 1/2 is that expected for the inner portion of 
NFW halos. It is tempting to conclude that halos are indeed NFW, and we will use this 
to compare the observed and predicted amplitudes £. An important caveat is that it is not 
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straightforward to match the radial scale of simulated halos to that of real galaxies. It is for 
this reason that the amplitude £ provides the most robust comparison: it uses all the data 
to set the velocity scale, in effect providing an integrated measure that is not sensitive to 
systematic effects in the way the innermost slope may be. 

The slope of the data is, if anything, too close to a constant 55 = 1/2. NFW profiles 
have this slope at small radii, though the exact slope the simulations predict is still debated 
(e.g., Ricotti 2003; Power et al. 2003; Navarro et al. 2004; Diemand et al. 2005). Irrespective 
of this detail, all simulations agree that the inner slope gradually transitions to a steeper fall 
off (becoming r~ 3 at large radii). The consequence is that there should be some curvature 
apparent in Vh(r) over the radial range sampled by the data. This is not obvious in the 
observations (Fig. [3]). This difference between simulations and data is sufficiently large that 
the difference between the various proposed analytic formulas is a rather small concern. The 
use of individually simulated halos (e.g., Hayashi et al. 2004; cf. de Blok 2005) in preference 
to the analytic approximations should not make a difference to the radii of interest here, 
where the deviations between the various realizations of the simulation results are small. 

In addition to the data, the upper panel of Fig. [3] shows the predicted shape and am- 
plitude of NFW rotation curves for a range of V^oo- Presuming galaxy halos span a range of 
mass comparable to that suggested by their luminosities, one would not only expect more 
curvature than observed, but also a wider spread in VJ, (compare, for example, the halos with 
V200 = 50 and 200 kms" 1 ). That the data concentrate in this plane more than anticipated 
might suggest that we may err in presuming that halo masses span the same range of mass 
as the observed luminosity, or that the simulated 0V200 relation built into the lines in the 
upper panel of Fig. [3] may not be precisely correct. 

One must bear in mind, of course, that the data only extend over a finite range of 
radii, and do not generally show the inevitable maximum and decline in V^. Though it is 
difficult to assign a specific halo to a particular galaxy, taken together the data are suggestive 
of an envelope composed of halos of a range of masses. If we had complete knowledge of 
the potentials out to arbitrary radii, we might see individual halos peeling away from this 
envelope as illustrated by the lines in the bottom panel of Fig. [3J 

For comparison with the data, we can ask what we would expect to find from a ACDM 
model observed in the same way as the data. We sample a series of halos with 50 < V200 < 300 
(e.g., the lines in the top panel of Fig. EJ). Model data points are allowed to fill the space 
available to them within the observed range of radii. The model points are uniformly sampled 
and weighted since we have no knowledge of the underlying halo mass distribution, other 
than that it gives rise to the broad observed range of luminosities and rotation velocities. 
Since no turn-over (V max of the halo) is observed in the data, the maximum radius to which 
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halos are sampled is truncated along the lower edge of the observed envelope, as defined by 
a line with slope 1/2 and intercept 1.33. The result of fitting the model halos is £ = 1.66 
and 33 = 0.48 for the vanilla ACDM parameters of Tegmark et al. (2004). For the third 
year WMAP parameters (Spergel et al. 2006), € decreases to 1.61. This normalization is 
too large, being consistent with the observed value only in the limit of zero stellar mass 
(£ = 1.62). 

The 0V200 relation can be adjusted in order to match the NFW model to the observed 
amplitude of the data. This is illustrated in the lower panel of Fig. [3] by lines representing 
NFW halos with parameters chosen to match the observed normalization <£. In effect, this 
provides an empirical estimate of the 0V200 relation. We can do this for the value of £ 
implied by any choice of T*. To get the zero point amplitude (£) for NFW, let us rewrite 
equation (jSJ) in log-log space. In the limit^l i^Owe have: 

logy, = \ [logy 2 oo + log[<?(c)] + log(V2)] + \ logr, (6) 

where 

9{C) = ln(l + c)-c/(l + c) (7) 

and the factor h = Ho/100 km s -1 Mpc -1 appears because V200 — ^-^200 when V200 is in 
kms" 1 and -R200 is i n kpc (NFW). Written this way, we recognize log 14 = £ + 03 logr with 
23 = 1/2. We can thus relate NFW parameters to the data through the observed £: 

log V 200 = 2<C - hg[g(c)} - log(/i/2) (8) 

(the data assume h = 0.75). In this fashion, we can generate the range of (c, V200) that are 
acceptable to the data (Fig. [TU|1 . Note that we do not specify the particular values of these 
parameters for any given galaxy. Rather, we constrain the range of C-V200 parameter space 
from which the observed halos may be drawn. 

The parameters of NFW halos in simulations are not arbitrary. There is a correlation 
between c and V200 (NFW) with fairly modest§| scatter (Jing 2000; Bullock et al. 2001; Colin 



2 The limit x — > might seem like a dubious approximation since we have excluded the data at r < 1 kpc. 
However, it is in fact quite a good approximation both theoretically (since x — r/i?2oo and we are concerned 
with data at a small fraction of the virial radius) and empirically (since the observed slope is close to limiting 
case of 1/2). 

3 Note that the scatter found by Colin et al. (2004) is half that found by Bullock et al. (2001). We adopt 
the more conservative (larger) scatter here. Adopting the smaller scatter would simply make the predicted 
bands in Fig. [TU] 2cr wide. 
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et al. 2004). Moreover, the typical concentration depends on the cosmology (NFW). For a 
10 12 M Q halo, McGaugh et al. (2003) found that the mean 

c = 1.88 + 23.9a 8 X .6 (9) 

where 

Xo.e = n^he^-^^ - 0.32(n; 1 - 1). (10) 

NFW found that the density of halos depends on the density of the universe at the time of 
collapse. This depends on Q m and the amplitude of the power spectrum at the appropriate 
scale, as encapsulated by the above formula. The concentration varies slowly with scale V200 
(NFW). 

For illustration, we adopt ACDM parameters from the Tegmark et al. (2004) and Spergel 
et al. (2006). We use the 6-parameter 'vanilla' fit of Tegmark et al. (2004), and the results of 
the 3 year WMAP-only fit from Spergel et al. (2006). While superficially these may appear 
similar, the concentrations of dark matter halos are quite sensitive to cosmic parameters, so 
there is a noticeable difference. 

The observed and predicted 0V200 bands are illustrated in Fig. [TUJ For the scatter in 
ACDM concentrations, we adopt a\ ogc = 0.14 (Bullock et al. 2001). Similarly, we illustrate 
the ±1 a range of observed halos (i.e., the band that contains 68% of the halos from the 
histograms in Figs. [6] and [7]). The predicted and observed bands do not match up well, 
having only a modest range of overlap at surprisingly small V 2 oo- 



5. Discussion 

The cusp-core problem in dwarf and LSB galaxies is well known. Inspection of Fig. [10] 
reveals a more general problem with the densities of the dark matter halos of all spiral galax- 
ies. The data are consistent with a swath of 0V200 space that does not parallel the predictions 
of ACDM. The acceptable regions of halo parameter space diverge as mass increases. 

This result can not be attributed to inadequate resolution or slit positioning errors. All 
of the data are derived from resolved two dimensional velocity fields. We have intentionally 
excluded data from the innermost region where the resolution might be questioned. Signif- 
icant non-circular motions could be seen directly in the velocity fields. The result does not 
change if we exclude^ even more of the inner data (see Fig. [Q imagine truncating Fig. [3] at 



4 This result is slightly less restrictive than the maximum halo acceleration limit (Brada & Milgrom 1999; 
McGaugh 2004) because we ignore here r = where NFW halos provide their maximum acceleration. 
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2 or 3 kpc). It is the velocity attributable to dark matter at well resolved intermediate radii 
that is lower than expected in ACDM. 

Let us examine the implications of this result, and possible ways out, starting with the 
theoretical expectations. Two fiducial ACDM models are illustrated in Fig. [TUJ The model 
of Tegmark et al. (2004) is a good representation of standard vanilla ACDM, predicting very 
nearly the same swath of 0-V200 parameter space as other plausible choices for cosmo logical 
parameters, going back to the early work of NFW. The WMAP three year results (Spergel 
et al. 2006) noticeably lower the predicted halo concentrations. This is due to slight changes 
in all of the relevant parameters. The dominant effect is from the lower power spectrum 
normalization a 8 = 0.90 — > 0.74, a result anticipated both by cluster (e.g., Rosati, Borgani, 
& Norman, 2002; Diego et al. 2003) and galaxy work (Zentner & Bullock 2002; McGaugh 
et al. 2003; van den Bosch, Mo, & Yang 2003). The tilt (n = 0.99 -> 0.95) and the lower 
matter density (fi m = 0.27 — ■> 0.24) also contribute. As we add in large scale constraints 
beside WMAP, the parameters tend_l to drift back towards those of Tegmark et al. (2004). 

The lower concentrations implied by the more recent WMAP results do help. Ignoring 
the trend with V200 for the moment, equation ([9]) implies c = 7.6 for the cosmic parameters of 
Tegmark et al. (2004), whereas the new WMAP data (Spergel et al. 2006) imply c = 5.9. This 
difference is quite pronounced in rotation curves. The revised WMAP results at least bring 
cosmic data into the right ballpark: McGaugh et al. (2003) found c ~ 5.7 for dark matter 
dominated galaxies, though it should be noted that this typical value ignores problematic 
cases with c < 2 or for which NFW simply does not fit. 

The problem here is not so much the mean concentration, but the run of concen- 
tration with halo mass. Concentrations are acceptable for low mass halos (those with 
^200 i$ 100km s -1 ). They become increasingly discrepant with increasing mass. By the 
time we reach V200 — 163 kms -1 (10 12 A / t ), the data imply c = 3. This is nominally the 
mass scale of the Milky Way (e.g., Wilkinson & Evans 1999; Sakamoto, Chiba, & Beers 
2003), and the low inferred concentration is consistent with both Milky Way constraints 
(Binney & Evans 2001; Flynn et al. 2006) and those for comparable galaxies (e.g, Weiner et 
al. 2001). To obtain such a low concentration from initial conditions would require absurd 
changes to cosmology: <7 8 ~ 0.3 or fi m « 0.05. 

This basic result is virtually independent of the choice of mass-to-light ratio estimator. 
The distribution of Vh is not as sharply peaked for V = 1 as it is for Q = 1, but the mean 
is nearly the same, so the same result is obtained. Worse, changing V (or Q or T) simply 
shifts the bands left (for V > 1) or right (for V < 1) in Fig. [TUJ It does little to alter the 



5 See |http://lam bda.gsfc.nasa.gov/product/map/current/parameters.cfm 
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differing slope of the observed and predicted 0V200 bands. 

In order to fix the mismatched slopes with T*, we would have to make the stellar mass 
indicator a systematic function of V^oo- Specifically, Y* would have to get progressively 
smaller as halo mass increases. This is quite contrived. There is no evidence to suggest that 
the stellar populations of massive spirals are systematically lighter than those of dwarfs. 
Indeed, such a contention can likely already be excluded by the need for a substantial stellar 
contribution to the mass in at least some high surface brightness galaxies (e.g., Debattista & 
Sellwood 2000; Weiner et al. 2001; Kassin et al. 2006; Battaglia et al. 2006). In contradiction 
to this, Jimenez, Verde, & Oh (2003) argue that rotation curve fits tolerate the expected 
concentrations for NFW halos. However, this comes at the expense of unrealistic T*. The 
/-band Y* distribution of Jimenez et al. (2003) is bimodal (their Fig. 7), with roughly 
half of the galaxies centered on Y* m 0. 2A4 Q /L Q . So they in fact find the same result, 
that mass-to-light ratios need to be quite low in order to be reconciled with the expected 
NFW concentrations. A velocity-variable IMF might achieve this, but it would introduce an 
unreasonable amount of scatter into the baryonic Tully-Fisher relation (McGaugh 2005). 

The mismatch between the slopes of the observational and theoretical bands in Fig. [TD] 
is caused by the near universality of the empirical dark matter halo found in §3. The data 
can be described, albeit crudely, by a single simple form (equation According to NFW, 
this should not be the case — halos should spread out (e.g., top panel of Fig. [3]), not clump 
together as observed. In this case, we would expect to see more data at high Vh at small r. 
While there are choices of Y* that have some data in this region, the bulk of the data do 
not reside there (Fig. Hj). 

The slope of the observed bands in Fig. [10] is the result of adjusting the NFW halos 
concentrations to match the data with 03 = 1/2. It is worth checking how this changes if we 
allow the slope to be free. A T* choice that leads to a shallower slope and higher amplitude 
is needed to reconcile the data with the model. 

Leaving both the slope and the intercept free in the fit, we obtain € = 1.43 and 23 = 0.58 
for V — 1. This steeper slope goes in the wrong direction, as it would make the V — 1 band 
in Fig. [TU] cut somewhat more steeply across the predicted band. Using the two dozen 
galaxies for which we have i^-band data and assuming a constant Y* = 0.7 M.q/L q (the 
V = 1 value for the mean color), we obtain CC = 1.53 and 23 = 0.49, consistent with 
Q — 1. Increasing the stellar mass makes things worse, maximum disk being the limiting 
case ((£ = 1.17 and 23 = 0.64 for T = 1): there is too little velocity attributable to dark 
matter at small radii. Lowering the stellar mass within reasonable bounds does not help 
either. For V = 1/2, (L = 1.59 and 23 = 0.49. This a band with precisely the same slope 
as illustrated; the change in €. only shifts the band a bit to the right. The case of T = 0.4 
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(Bottema 1993) gives results indistinguishable from V = 1/2. To push the stellar mass 
low enough to get into the right ballpark, we need V < 0.3 (C = 1.65 and 03 = 0.44 for 
V = 0.3). This corresponds to an IMF that is a factor of 5 less massive than Salpeter. 
This is below the practical minimum in stellar mass as determined from direct integration 
of the observed IMF (Kroupa 2002). It is also grossly inconsistent with the coherence of the 
mass discrepancy-acceleration relation (McGaugh 2004) and baryonic Tully-Fisher relation 
(McGaugh 2005), the presence of baryonic features observed in rotation curves (Palunas & 
Williams 2000; Sancisi 2004; Kassin et al. 2006; Battaglia et al. 2006), observed velocity 
dispersions (Bottema 1997; Kregel, van der Kruit & Freeman 2005), and hydrodynamical 
flows (Weiner et al. 2001; Kranz, Slyz & Rix 2003). There appears to be no plausible choice 
of Y* that is consistent with cosmological expectations and NFW, at least for halos with 
V200 comparable to the observed circular velocity. 

We have so far implicitly assumed that halo mass scales linearly with baryonic mass. 
This need not be the case. One possibility to consider it that the observed luminosity function 
of galaxies descends from parent dark matter halos spanning a smaller range of mass. This 
would require a systematic variation in the fraction of baryons that form the luminous disk 
(McGaugh & de Blok 1998; Okamoto et al. 2005), but this may already be required to match 
the observed Tully-Fisher relation (van den Bosch 2000; Mo & Mao 2004; see also Navarro 
& Steinmetz 2000; Eke et al. 2001; Portinari & Sommer-Larsen 2006; Dutton et al. 2006). 

If we drop the assumption that total mass scales with baryonic mass, it is tempting 
to conclude that galaxies may be born from halos in the region allowed in Fig. [IUJ This 
would mean that even massive galaxies with observed rotation velocities in excess of 200 
kms" 1 (ranging all the way up to 300 kms -1 ) would have their origins in relatively modest 
halos. From Fig. [IOJ we need V200 < 100km s" 1 for vanilla ACDM, and V 2 qq < 180km s -1 
for WMAP third year parameters. Galaxies would presumably gain their larger observed 
velocities during the collapse of their massive disks (see below). 

Examination of Figs. [3] and H] reveals that halos with V200 ~ 100km s" 1 are about right 
to explain much of the data. Unfortunately, such low mass halos do not provide a complete 
explanation. For one thing, there need to be enough baryons available in the halo to form 
the observed galaxy. Assuming a universal baryon fraction of 0.17 (Spergel et al. 2006), and 
the most favorable case that all baryons within i?2oo collapse to form the visible galaxy, a 
halo with V200 — 100km s -1 has enough baryons to make a galaxy with Vf = 167km s _1 . 
This assumes Q = 1 mass-to-light ratios and the baryonic Tully-Fisher relation of McGaugh 
(2005). There is some room to play with the choice of Y*, but this does little more than 
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assuage the extreme assumption that all baryongj collapse to form the observed galaxy. It is 
therefore rather difficult to ascribe the known spiral galaxies with 200 < Vf < 300km s -1 to 
NFW halos of low mass but normal concentration, at least for vanilla ACDM. This constraint 
is relaxed for WMAP 3 year parameters, but even here a halo with V200 = 180km s -1 can 
give at most enough baryons for a galaxy with Vf = 260km s _1 (again taking Q — 1). 

Looking in detail, NFW halos obeying the expected C-V200 relation have more curvature 
than observed, predicting too high a velocity at small r and too small a velocity at large 
r relative to the data. Even if we overlook the over-prediction at r ~ 5 kpc, the under- 
prediction at r > 20 kpc is quite serious. The compressive formation of a massive disk can 
raise Vh, but this process is much more effective at small radii than at large radii. In order 
to explain the high observed dark halo velocity at large radii, it seems necessary to have 
initial halos with at least V200 ~ 200km s -1 . Such halos grossly over-predict Vh at small r if 
they obey the vanilla 0V200 relation. It may be tempting to simply ignore the data at large 
radii, as there are relatively few with r > 30 kpc. However, there do exist galaxies outside 
this sample (e.g., Noordermeer et al. 2005; Spekkens & Giovanelli 2006) with outer rotation 
velocities approaching 300 kms -1 . These galaxies show no sign of a turn down in their outer 
velocities, and obey the same scaling relations as those in the present sample, so there is 
no reason to expect them to be any different. The need for large Vh at large radii is thus a 
serious problem. 

So far, we have compared the data to the ACDM predictions for primordial dark matter 
halos. It is possible that the process of galaxy formation alters the halos. If so, the difference 
between the predictions and observations might be interpreted to show what changes occur 
in this process. One should bear in mind that cuspy mass distributions seem to be a generic 
result of cold collapse (e.g., Huss et al. 1999), and that, once established, it is very hard to 
unbind the mass in the deep potential wells of cuspy halos. Achieving the desired effect is 
thus a tall order from the outset (Sommer-Larsen & Dolgov 2001; Eke et al. 2001). 

The primary physical effect which should occur during galaxy formation is the adiabatic 
contraction of the dark matter halo in response to the collapse of the baryonic disk (e.g., 
Blumenthal et al. 1986; Barnes & White 1984; Gnedin et al. 2004; Sellwood & McGaugh 
2005; Choi et al. 2006). The effect of adiabatic contraction is to further concentrate the dark 
matter. This goes in the wrong direction to help with the problems discussed here. However, 
there are two distinct aspects of adiabatic contraction to consider. The slope 7 will steepen 
from its initial value; this goes in the wrong direction. The velocity provided by the dark 
matter also increases. While this might also seem to act in the wrong direction, it means 



6 In ACDM, the virial radius is closer to i?ioo than -R200, further reducing the number of available baryons. 
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that a galaxy of a given rotation velocity may form in a halo with lower primordial V^oo- As 
discussed above, this can help, at least in principle. Unfortunately, it seems unlikely that 
the effect can be large enough to help at large radii without simultaneously doing far too 
much at small radii. Disks heavy enough to have a substantial impact on outer velocities 
grossly over-predict inner velocities if the starting point is an NFW halo (Abadi et al. 2003; 
Sellwood & McGaugh 2005). To make matters worse, Dutton et al. (2006) argue that a 
large difference between the observed and initial halo velocity scales would make reconciling 
the normalization problem between Tully-Fisher and the luminosity function practically 
impossible. Low mass halos might also be inconsistent with constraints from weak lensing 
(Smith et al. 2001; Mandelbaum et al. 2006). 

After adiabatic contraction, there may be some subsequent process that further alters 
the halo. Perhaps the most common example is feedback from star formation. Quite generi- 
cally, the effects of feedback are believed to be stronger in lower mass galaxies owing to their 
shallower potential wells (e.g., Thacker & Couchman 2000). Such an effect is unlikely to be 
sufficiently strong (Mac Low & Ferrara 1999; Gnedin & Zhao 2002) and acts preferentially 
on the wrong end of the mass function. It is the concentrations of high mass galaxies (if 
not necessarily the slope 7) that are most deviant from the predictions of ACDM. If the 
solution to this problem lies in the processes of galaxy formation, it likely resides in a rather 
different direction (see, for example, El-Zant et al. 2001; Weinberg & Katz 2002; Merritt 
et al. 2004; Ma & Boylan-Kolchin 2004; Mashchenko, Couchman, & Wadsley 2006; Tonini, 
Lapi, & Salucci 2006; see also Sellwood 2006). 

Another approach would be to alter the initial conditions. A turndown in the power 
spectrum on galaxy scales can reduce halo concentrations. However, such a solution would 
require a rather strange dip in the power spectrum. Concentrations at small scales are not 
unreasonable, and also make sense on larger (cluster) scales (e.g., Vikhlinin et al. 2006). It 
is the intermediate scale of massive spirals, giant ellipticals (Romanowsky et al. 2003; Pierce 
et al. 2006), and groups (Karachentsev 2005) where the halo concentration prediction^] of 
ACDM simulations fares worst. 

Another possibility is to modify the nature of dark matter. The cusp-core problem has 
already motivated consideration of various possibilities, such as warm (Bode, Ostriker, & 
Turok 2001), fuzzy (Hu, Barkana, & Gruzinov 2000), self-interacting (Spergel & Steinhardt 
2000), and meta-cold (Strigari et al. 2006) dark matter, among others (e.g., Piazza & Mari- 



7 The prediction of a steep inner slope, 7 « 1 — the shape rather than the amplitude — continues to 
suffer most clearly in LSB systems (Simon et al. 2005; Kuzio de Naray et al. 2006; Goerdt et al. 2006; 
Sanchez- Salcedo et al. 2006). 
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noni 2003; Lee & Lee 2004; Blanchet 2006). It is not immediately clear how any of these 
proposals fare in the current context. However, many proposals seem to be meant to produce 
dark matter halos with a core. This in itself is not a complete solution. One must get the 
density right, and do so on all relevant scales. 

6. Conclusions 

We have examined the rotation curves of a large sample of spiral galaxies in order to 
constrain the velocity that can be attributed to dark matter. Taken together, the dark halo 
velocity can be approximated by 

logVh = £ + <Blogr 

with 23 ~ 1/2 and £ = 1.47l2;J|. The value of <£ depends on the choice of stellar mass-to-light 
ratio, with the uncertainty corresponding to the full range of possibilities from minimum to 
maximum disk. While the dark matter halos of individual galaxies must differ, this relation 
gives a tolerable approximation over the observed range. 

In the context of ACDM, the densities of halos depends on cosmology. For reasonable 
cosmological parameters, we expect £ > 1.6. This is only consistent with the data in the 
limit Y* — > 0, which is excluded by a wide variety of independent lines of evidence. The 
problem is not just that rotation curve data prefer halo models with a core over those with 
a cusp. The amplitude of the velocity provided by dark matter is too low at all observed 
radii. 

We have emphasized here the importance of spiral galaxy dark matter halo densities 
being substantially lower than anticipated by ACDM. For halos of 10 12 A-f©, the density dis- 
crepancy is nearly an order of magnitude (see also Sellwood 2006). Getting the halo density 
in the right ballpark would be an important step towards a successful theory. Ultimately, 
one must explain in detail the observed coupling between dark and baryonic components 
(e.g., Sancisi 2004; McGaugh 2004). 
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Fig. 1. — The radial distribution of the data from all 60 galaxies. Data with r < 1 kpc (un- 
shaded portion of histogram) are excluded from further analysis. The dotted lines illustrate 
the range 4 < r < 22 kpc containing 68% of the data, with equal numbers on either side of 
the median (9 kpc). These data at intermediate radii dominate the signal and subsequent 
result. 
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Fig. 2. — The rotation curve and mass model of NGC 2403 in log-log space. The Hi data 
are from Begeman (1987; see also Fraternali et al. 2002), Ha (i) from Blais-Ouellette et al. 
(2004), and Ha (ii) from Daigle et al. (2006). The Ha data are shown to illustrate the degree 
to which independent data sets agree, but only the HI data are used in the analysis. The 
solid lines show the mass model Vb of the baryons for two choices of Y*, T = 0.9 (upper) and 
V = 1 (lower). The dashed lines show the resulting halo rotation curves Vh for each Y*. The 
upper dashed line corresponds to the lower baryonic line and vice-versa. The dotted line at 
r = 1 kpc demarcates where the innermost data are ignored. 



ACDM c-V 2QQ 



10° 10 1 10 2 



R (kpc) 




10° 10 1 10 2 

R (kpc) 

Fig. 3. — The rotation velocity Vh attributable to the dark matter halos of 60 galaxies 
(points) for Q = 1 stellar masses. Taken together, these data form a narrow sequence that 
is well described by logV/j = 1.47 + 0.5 log R (straight line). Also shown are galaxy mass 
NFW halos with V200 — 50, 100, 200, and 300 kms -1 (curved lines). In the top panel, these 
are drawn for the concentrations predicted for the vanilla ACDM parameters of Tegmark et 
al. (2004). In the lower panel, the c — V200 relation has been adjusted to match the data. 
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Fig. 4. — Dark halo velocities, as per Fig. [3l The data are represented here as a grayscale to 
show the density of the data. Darker shades representing regions of higher density; lighter 
shades those of low density. In regions of very low data density, individual points are shown 
as dots. While the eye may be attracted to outliers in Fig. [3], here it is clear how sharply 
defined Vh (r) is. Several choices of mass-to-light ratio are illustrated: maximum disk (r = 1, 
top left); stellar population synthesis (V = 1, top right); MOND (Q — 1, bottom left - the 
same as in Fig. [3]); and sub- maximal (Bottema 1993) disks (r = 0.4, bottom right). 
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Fig. 5. — The distribution about the line fit to the Q = 1 case for maximum disk mass-to- 
light ratios (T = 1). In essence, this is the cross-section of the corresponding pane in Fig. HI 
All the data are shown in the first panel, with subsequent panels showing the data falling 
in radial bins containing approximately equal numbers of data points. Scatter causes some 
individual points to imply a negative dark mass. These are placed in the bin at the left edge 
of the plot. The variation of the mean SlogV illustrates how the characteristic amplitude of 
the dark matter velocity tt shifts with T*, becoming lower in this case. The radially binned 
plots illustrate how the distribution varies. That the centroid in these bins does not vary 
much shows the approximate constancy of the slope 23. 
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Fig. 6. — As per Fig. [51 but for population synthesis stellar mass estimates (P = 1). 
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Fig. 7.— As per Fig. Eg but for MOND (Q 
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Fig. 8. — As per Fig. [51 but for sub-maximal disks with F = 0.4. 
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Fig. 9. — Variation of the scatter (left) and mean zero point velocity (AtC from Q — 1; right) 
with stellar mass prescription scaling factor (Q: solid line; V: dotted line; V: dashed line). 
A clear minimum in the scatter occurs at Q — 1. The variation for Q > 1 is not smooth as 
individual galaxies are excluded as the prescription causes them to exceed maximum disk. 
Note that for low stellar masses (Q, V < 0.8, T < 0.4) the different prescriptions give 
indistinguishable mean (£. 
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Fig. 10. — The dark matter halo concentration-virial velocity plane. Blue (V = 1) and gray 
(Q = 1) bands show the region consistent with the data for those choices of mass-to-light 
ratio estimator. Orange and red bands show the region that NFW halos are predicted to 
occupy for the cosmological parameters of Tegmark et al. (2004) and Spergel et al. (2006), 
respectively. In both cases, the width of the bands is ±lo~. The predicted bands assume an 
intrinsic scatter of a c = 0.14 (Bullock et al. 2001). The width of the predicted bands would 
shrink and the region where they overlap would be reduced if we adopted the smaller scatter 
(ct c = 0.11 for relaxed halos) advocated by Colin et al. (2004). The predicted concentrations 
overlap with the data over a rather narrow range. This range would of course shrink for a 
smaller intrinsic scatter in c. Theory and data diverge with increasing V^oo- The divergence 
is most pronounced for the range of halo masses typically assumed to be associated with L* 
galaxies. 



